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Abstract 

A new type of the single-molecule magnet [Mn403Cl4(02CEt)3(py)3] forms dimers. 
Recent magnetic hysteresis measurements on this single-molecular magnet revealed 
interesting phenomena: an absence of quantum tunneling at zero magnetic field and 
tunneling before magnetic field reversal. This is attributed to a significant antifer- 
romagnetic exchange interaction between different monomers. To investigate this 
system, we calculate the electronic structure, magnetic properties, intramolecular 
and intermolecular exchange interactions using density-functional theory within the 
generalized-gradient approximation. Our calculations agree with experiment. We 
also calculate vibrational infrared absorption and Raman scattering intensities for 
the monomer which can be tested experimentally. 
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1 Introduction 



Single-molecule magnets (SMMs) are three-dimensional arrays of identical 
nanoscale molecules, each of which consists of several transition metal ions 
surrounded by organic ligands. These SMMs are a promising building block for 
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magnetic storage devices and for investigating macroscopic quantum tunneling[l] 
and quantum decoherence.[2] A prototype of these SMMs is [Mni20i2 
(CH3COO)i6(H20)4]-2(CH3COOH)-4(H20) (hereafter Mni2)[3] which has an 
effective ground-state spin of 5" = 10. Magnetic hysteresis measurements 
on the SMM Mni2 showed quantum tunnehng between spin-up states and 
spin-down states. [4,5,6,7] This has been successfully explained by magnetic 
anisotropy parameters only. However, recent measurements on a new type of 
SMM, [Mn403Cl4(02CEt)3(py)3] (hereafter Mn4)[8,9] where Et=CH2CH3 and 
py=NC5H5=pyridine, showed qualitatively different tunneling behavior [10]: 
quantum tunneling occurred before the magnetic field was reversed and no 
quantum tunneling was observed at zero field, which is in contrast to other 
SMMs like Mni2. It was speculated that this behavior was caused by sub- 
stantial antiferromagnctic exchange interactions between different monomers 
within a dimer.[10] Suppose that magnetic moments of all molecules of the 
SMM Mn4 are aligned due to strong external magnetic field. If the direction 
of the external magnetic field is reversed adiabatically through zero field, then 
some molecules may change their magnetic moments before the external field 
reversal. This is because the antiferromagnctic exchange interactions facilitate 
monomeric magnetic moment flips prior to the external fleld reversal. This in- 
teraction also prevents two monomers from simultaneously flipping their mag- 
netic moments at zero field. In contrast, it is widely accepted that for Mni2 
there are no appreciable exchange interactions between different monomers. In 
this sense, further investigation of the SMM Mn4 is worthwhile. Additionally, 
the SMM Mn4 consists of the same kinds of Mn spins as those for Mni2 but 
it is much simpler than Mni2 since the number of Mn spins is smaller. There- 
fore the SMM Mn4 can be used to obtain some insights into issues remaining 
in the SMM Mni2 such as clarifying the origins of symmetry-breaking terms 
in the effective single-spin Hamiltonian[ll,12] and possible spin-orbit-vibron 
interactions. [13] 

A dimeric form of the SMM Mn4 is obtained by inversion symmetry of the 
threefold symmetric monomer shown in Fig. 1. Since dimers are well sepa- 
rated from each other, interactions between different dimers are negligible so 
that they are not considered in our calculations. The magnetic core of the 
Mn4 monomer comprises three ferromagnetically coupled Mn^"*" (S = 2) ions 
coupled antiferromagnetically to the remaining Mn^+ {S = 3/2) ion leading 
to a total ground-state spin of S = 2 x 3 — 3/2 x 1 = 9/2. The core has a 
similar cubane structure to the inner core of the SMM Mni2, although in Mni2 
there are four Mn'^^ ions. In this work, we investigate the electronic structure, 
magnetic properties, and vibrational modes of the SMM Mn4 using density- 
functional theory (DFT).[14] In Sec. 2.1, we present calculations of the Mn4 
monomer such as optimized geometries, intramolecular exchange interactions 
between Mn spins, and magnetic anisotropy barriers for the ground-state and 
excited-state manifolds. In Sec. 2.2, we discuss binding energy and exchange in- 
teractions between different monomers within a dimer. In Sec. 2.3, we present 
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calculated infrared and Raman intensities for vibrational normal modes. In 
Sec. 3, we present our conclusion. 



2 DFT calculations and Discussion 

We have used all-electron Gaussian-orbital-based Naval Research Laboratory 
Molecular Orbital Library (NRLMOL) [15,16,17,18,19,20,21] within the Perdew- 
Burke-Ernzerhof (PBE) generalized gradient approximation (GGA)[22] along 
with full basis sets for all atoms and fine mesh. [23] In our analysis, we have 
considered monomers and dimers that are terminated by both H and by the 
CH2CH3 radicals found in the experiments. We have constructed a dimer com- 
prising two conformers of the monomeric units. Our study showed that both 
conformcrs are stable. Our vibrational analysis on a monomer terminated by 
H confirms that there are no unstable modes. The two conformers have slightly 
different arrangements of the pyridine ligands. The first conformer was found 
by the density-functional-based optimization of the hydrogenated geometry. It 
is called the computationally determined conformer (CDC). The second con- 
former was obtained by the experimental x-ray data with hydrogen positions 
corrected. This is called the experimentally determined conformer (EDC). 



2.1 Monomeric electronic and magnetic properties 

Details of finding the CDC monomer were discussed clscwhcrc[24] so here they 
are briefiy summarized. The Mn4 monomer has threefold symmetry about the 
direction of the central bridging Cl-Cl bond in the dimer geometry shown in 
Fig. 1 (or the direction of connecting Mn^"*" and CI in the cubane). To find an 
initial geometry, we optimize a pyridine ring and a cubane separately. Then 
the initial geometry for the monomer is relaxed using NRLMOL with the CI 
atom fixed to reproduce the experimental Cl-Cl distance upon dimerization. 
Relaxation continues until maximum forces between atoms are no larger than 
~ 0.001 hartree/bohr. A geometry for the dimer can be found by inversion 
symmetry of the optimized monomer geometry. Charges and magnetic mo- 
ments for the Mn atoms from the CDC monomer agree well with those from 
the EDC monomer. The total magnetic moment is Q/x^, which agrees with ex- 
periment. The energy gap between the minority lowest unoccupied molecular 
orbital (LUMO) and the majority highest occupied molecular orbital (HOMO) 
is much larger than the thermal energy shown in Table 1, so the total mag- 
netic moment of the ground state is stable. Calculated electronic density of 
states (not shown) confirms that the three Mn^+ spins are antiferromagneti- 
cally coupled to a Mn^+ spin. 
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Since both the CDC and EDC monomers are stable, we choose one of them, for 
example, the optimized CDC monomer to examine intramolecular exchange 
interactions. Considering the symmetry of the monomer and the fact that 
there are two kinds of Mn spins, we recognize that there are two types of ex- 
change interactions between Mn spins within a monomer as shown in Fig. 2. 
To calculate the exchange coupling constants, we assume that magnetic mo- 
ments of all Mn spins are aligned along a particular direction and of Ising 
type. So a Mn^"*" spin can have either = +2 or = —2 and a Mn''"'" 
spin can have = +3/2 or = —3/2. Thus we can construct three dif- 
ferent spin configurations other than the ground-state shown in Table 2. For 
example, Ms = 15/2 {Mg is an eigenvalue of the z component of the total 
spin operator S) is achieved by flipping a Mn^'^ spin from the ground-state 
5" = 9/2: = 2 2 2 3/2 = 15/2. Notice that all of the examined spin 
configurations are not eigenstates of S"^. We set up geometries of the three spin 
configurations using the optimized ground-state geometry with corresponding 
Ising-type spin arrangements and proper magnetic moments. Then we mini- 
mize self-consistently the energies of the spin configurations. Comparing the 
energies of the three configurations (refer to Table 2) with the energy gap 
between the minority LUMO and the majority HOMO, we find that the ex- 
amined spin configurations have an order of magnitude smaller energy than 
other kinds of spin excitations such as moving one majority spin to the unoc- 
cupied minority orbital. So we call those spin configurations low-energy spin 
excitations. Since there are three unknowns (the background energy, Eq, the 
two types of exchange constants, Ji and J2) and four equations to solve, we 
can calculate the exchange constants by the least-square-fit (LSF) method. 
Our calculated values of Ji and J2 are presented in Table 3. 

Differences between the calculated energies and the LSF-determined values 
for the ground state and low-energy spin excitations are quite small as shown 
in Table 2. The exchange constant between Mn^"*" spins, Ji, is confirmed to be 
ferromagnetic, and the exchange constant between Mn^"*" spins and a Mn'^+ 
spin, J2, is antiferromagnetic but larger than Ji. From exact diagonalization of 
Heisenberg exchange Hamiltonian, we confirm that the ground-state manifold 
has S = 9/2, that the first excited-state manifolds are doubly degenerate with 
S — 7/2, and that the second excited-state manifold has -S" = 11/2. The energy 
gaps between the ground-state manifold and the excited-state manifolds are 
given in Table 3. The large energy gap between the ground-state and the 
first-excited manifold also supports that the ground-state S = 9/2 manifold is 
stable. To check how good the calculated values of Ji and J2 are compared to 
experiment, we calculate the effective moment per molecule, /leS, at magnetic 
field of 1 T as a function of temperature T as follows: 
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where x is susceptibility, /c^ is the Bohzmann constant, and jiB is the Bohr 
magneton. If we use the calculated vahies of Ji and J2 given in DFT(l) of 
Table 3, then we obtain overestimated effective moments (dashed curve in 
Fig. 3) at high temperatures compared to the experimental data. If we use 
half-reduced values of Ji and J2 from the DFT-calculated ones, then the cal- 
culated effective moment (solid curve in Fig. 3) agrees well with experiment. 
Notice that the reduced values of Ji and J2 are quite close to the experi- 
mentally extracted values in Table 3. This supports the argument that DFT 
calculations often overestimate the exchange interactions between atoms be- 
cause of imperfect treatment of self-interaction of the Coulomb potential. 

We calculate the magnetic anisotropy barrier (MAE) in zero magnetic field for 
both the CDC and EDC monomers for the ground state S = 9/2, assuming 
that the spin-orbit coupling mainly contributes to the barrier. Following the 
procedure explained in Ref . [25] , we find that the Mn4 monomer has the easy 
axis along the threefold axis, in agreement with experiment. The calculated 
anisotropy barrier is about 11.3 K for the CDC monomer, 11.6 K for the EDC 
monomer, and 10.9 K for the hydrogenated EDC monomer. All of these values 
are close to the experimentally measured value of 14.4 K. To investigate the 
contribution of each Mn spin to the barrier, we project the total barrier onto 
each Mn spin and calculate its projected anisotropy. The three Mn^"*" spins 
have the easy axes along the x, y, and z axis respectively. The Mn^"*" spin 
has the easy axis along the (111) direction. Most of the anisotropy originates 
from the Mn^+ spins due to Jahn-Teller distortion, while the contribution of 
the Mn^+ spin is insignificant. The same tendency occurs for the SMM Mni2. 
We also calculate the magnetic anisotropy barriers for the low-energy spin 
excitations shown in Table 2. The barriers are essentially the same as that for 
the ground state S = 9/2. 



2.2 Dimeric properties 

For a dimeric form, the distance between the central bridging Cl-Cl bond is 
kept as the experimental value of 3.86 A unless specified. We calculate the 
binding energy of the dimer by subtracting the dimer energy from twice the 
monomer energy. We find that the dimer is stable for both the CDC and 
EDC. For the CDC (hydrogenated EDC) dimer, the binding energy is about 
0.16 eV (0.45 eV).[24] The magnitude of the binding energy suggests attractive 
electrostatic interactions between different monomers. We find that the electric 
dipole moment for each CDC (hydrogenated EDC) monomer is as large as 2.28 
(1.91) in atomic units and that it points towards the CI atom in the cubane 
along the threefold axis. Since the dimer has inversion symmetry, the electric 
dipole moment for the dimer vanishes. The discrepancy between the binding 
energy for the CDC and that for the EDC arises from our substitution of 
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CH2CH3 for H and slightly different geometries of pyridine rings for both 
conformers. 

To calculate the exchange coupling constant Jmter between different monomers 
within a dimer, we make assumptions that a monomer is an ideal S — 9/2 
object and that its effective spin is aligned along the easy axis and of Ising 
type (either = +9/2 or = —9/2). Then we calculate sclf-consistently 
the energies of ferromagnetic (parallel monomeric spins) and antiferromag- 
netic (antiparallel monomeric spins) configurations of the dimer, and take the 
difference S between the two energies. We find that the antiferromagnetic con- 
figuration is favored, and that from S = 2 Jintcr(9/2)^ the exchange constant 
Jintcr is about 0.24 K for the CDC and 0.27 K for the EDC.[24] This can be 
compared to the experimentally measured value of 0.1 K. From the analysis 
of the exchange constants within a monomer, this overestimated value of Jinter 
is not so surprising. We also find that Jinter is quite sensitive to the central 
bridging Cl-Cl bond length. Our calculations show that Jinter increases expo- 
nentially with decreasing the Cl-Cl distance. More detailed analysis was given 
in Ref.[24]. Since the experimental results cannot distinguish exchange inter- 
action from magnetic dipole-dipole interaction, we estimate the contribution 
of the magnetic dipole-dipole interaction to the total interactions between dif- 
ferent monomers in a dimeric form. Since the dipole moment is along the easy 
axis and the two monomers are located along the easy axis, the ferromagnetic 
configuration is favored. But the difference in the dipole-dipole interaction be- 
tween the ferromagnetic and antiferromagnetic configuration is about 0.28 K, 
which is neghgible compared to the exchange energy difference of 10 K. This 
was calculated considering two large spins of S = 9/2. Including the internal 
structure of a monomer does not substantially change the energy difference in 
the dipole-dipole interaction. 



2.3 Vibrational analysis 

We calculate vibrational normal modes for the EDC monomer with CH2CH3 
substituted by H within the harmonic oscillator approximation. There are a to- 
tal of 168 normal modes (= 56x3), six of which are translational and rotational 
modes. All of these modes are vibrationally stable, which is also corroborated 
by our frozen phonon calculations. Based on the normal coordinates obtained 
from the vibrational analysis, we investigate which modes are infrared (IR) 
and/or Raman active. IR absorption intensity is proportional to the square of 
change of electric dipole moment with respect to displacement along the nor- 
mal coordinates. Raman scattering intensity is proportional to square of the 
change of polarizability with respect to displacement along the normal coordi- 
nates. We calculate the IR absorption and Raman scattering intensities follow- 
ing the method suggested in Ref. [21,26,27,28]. First, we discuss the IR spectra 
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shown in Fig. 4. Projection of the total IR intensities onto each component in 
the monomer, indicates which component contributes to the spectra at cer- 
tain modes. Overall, we find four strong IR intensities at 543 cm~^, 1316 cm~^, 
1547 cm"-*^, and 2942 cm~^. For low-frequency modes below 444 cm~^ all parts 
of the monomer contribute. At 543 cm~^, the dominant contributions are from 
and Mn. The contributions of Mn and are limited to low-frequency 
modes (below 602 cm^^). From 671 cm~^ the contribution from pyridine rings 
becomes important. At 1316 cm~^, 1547 cm~^, 1563 cm~^, and 2942 cm~^ the 
contributions are from O2CH and pyridine rings. For high-frequency modes 
the vibrations of O2CH and pyridine rings provide dominant IR intensities. 
Second, we discuss the Raman spectra shown in Fig. 5. The Raman intensi- 
ties are strong mostly at very low vibrational frequencies, in contrast to the 
IR intensities that are strong at high frequencies. The strongest Raman in- 
tensity is at 21 cm~^, and the next largest peak is at 35 cm~^. These modes 
are mainly due to vibrations of pyridine rings as well as small contributions 
from CI, Mn, and O2CH. At 76 cm"-*^ and 96 cm~^ the dominant contributions 
are from CI and smaller ones from pyridine, O2CH, and Mn. At 111 cm~^, 
168 cm~^, and 212 cm~^, the contributions are from all components except 
for 0^~. At 514 cm~^, 543 cm"^, and 602 cm~^ the contributions are from 
and Mn only. At 1019 cm^^, only pyridine rings contribute, and at 2942 cm~^ 
(not shown) O2CH and pyridine rings contribute. At 3115 cm~^, 3132 cm~^, 
and 3154 cm^^ (not shown) only pyridine rings contribute. Compared to the 
IR spectra, the contributions of CI atoms are noticeable for low-frequency 
modes and those of O^^ appear at between 514 cm^^ and 602 cm^^ only. For 
high-frequency modes, the contributions are from pyridine and O2CH only. 



3 Conclusion 



We have found optimized geometries for the Mn4 monomer and dimer us- 
ing DFT, and calculated exchange interactions between Mn spins within the 
monomer and magnetic anisotropy barriers for the ground state S = 9/2 and 
the low-energy spin excitations. We also have calculated the binding energy 
and the monomer-monomer exchange interaction and magnetic dipole-dipole 
interaction within a dimer. Most of our DFT calculations are in good agree- 
ment with experiment, although the calculated intermolecular/intramolecular 
exchange interactions are somewhat overestimated compared to the experi- 
mental values as is common for DFT calculations. Our vibrational analysis 
on the hydrogenated EDC monomer rules out a possibility of breaking the 
threefold symmetry by stable soft vibrational modes in the SMM Mn4. The 
calculated IR and Raman spectra can be checked by experiments. 
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Table 1 

Energy gaps between the highest occupied molecular orbital (HOMO) and the lowest 
unoccupied molecular orbital (LUMO) for majority and minority spins for the SMM 
Mn4 monomer^ 



maj HOMO - maj LUMO 


min HOMO - min LUMO 


1.02 eV 


2.42 eV 


min LUMO - maj HOMO 


maj LUMO - min HOMO 


1.17 eV 


2.28 eV 



Table 2 

Eigenvalues of the z component of the total spin operator, M^, spin configurations, 
their Ising exchange energies, calculated energies relative to the ground state energy, 
A£J, differences between the calculated and the least-square-fitted energies, and 
calculated magnetic anisotropy barriers (MAE) per monomer for the ground state 
and the low-energy spin excitations. For the spin configurations, the first three 
symbols are for the three Mn"^+ spins ("+" denotes = +2) and the fourth one 
is for the Mn^"*" spin ("-|-" denotes = +3/2). Here Eq is the background energy, 
Ji and J2 are the intramolecular exchange constants indicated in Fig. 2. These 
calculations were performed for a computationally determined conformer (CDC). 



Ms 


spin conf. 


Ising energy 


(eV) 




MAE/monomer (K) 


9/2 


+ + + - 


Eo + 12 Ji - 9 J2 





-2.7 X 10-4 


11.3 


1/2 


+ + 


Eq - 4Ji - 3 J2 


0.14 


8 X 10^4 


11.4 


7/2 


+ + - + 


Eq - 4Ji + 3 J2 


0.217 


-8 X 10""^ 


11.8 


15/2 


+ + + + 


Eq + 12 Ji + 9J2 


0.236 


2.7 X 10-4 


10.9 



Table 3 

Calculated intramolecular exchange constants, energy gap between the ground-state 
S = 9/2 and the first excited doubly degenerate 5" = 7/2 manifold, E{S = 7/2), 
and energy gap between the ground-state and the second-excited S = 11/2 mani- 
fold, E{S = 11/2), in comparison with the experimental results. [8,10] The negative 
sign means ferromagnetically coupled. DFT(l) denotes our calculated values, while 
DFT(2) denotes values obtained by reducing our calculated values by half. The last 
two columns are the binding energy and intermolecular exchange constant within a 



CDC dimer. 





Ji(K) 


J2(K) 


E{S = DiK) 


E{S = ^) (K) 


Binding energy (eV) 


•Winter (-^) 


DFT(l) 


-44 


152 


490 


836 


0.16 


0.24 


DFT(2) 


-22 


76 


245 


418 






Exp 


-25[8] 


60 [8] 


238 [8] 


330 [8] 




0.1[10] 
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Fig. 1. (Left) Cubane-like magnetic core of the Mn4 monomer. The four arrows 
represent Mn spins and O and CI atoms are labeled. The magnetic core consists of 
three ferromagnetically coupled Mn^"'" spins {S=2) coupled antiferromagnetically to 
one spin {S=3/2) ion leading to a total spin of 5 = 9/2. (Right) Mn4 dimer 

geometry. The dimer is formed by inversion of the threefold symmetric monomer. 
The distance between the two central CI atoms marked as the dotted line was 
measured to be 3.86 A. The threefold axis is along this Cl-Cl bond. 




Fig. 2. Schematic diagram of intramolecular exchange interactions within a 
monomer. Three Mn^+ (S = 2) spins are denoted as empty circles and one Mn^+ 
{S = 3/2) spin as a filled circle at the four vertices of the tetrahedra. Each spin is 
numerically labeled at each vertex. Ji is the exchange interaction between Mn'^+ 
spins and J2 is between Mn'^"^ and Mn'^+ spins. 
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132 198 
Temperature (K) 



Fig. 3. Calculated effective moment /XeS per monomer. The dashed, solid, and 

dot-dashed cTirvcs arc obtained using the calculated values of Ji and J2, a half 
of the calculated values, and a quarter of the calculated values respectively. The 
solid line is close to the experimental data. [8] 
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Fig. 4. Calculated total and projected infrared (IR) vibrational spectra for the 
experimentally determined conformer (EDC) with CH2CH3 replaced by H. From 
the top panel, shown are the projected IR active density of states onto Mn, CI, 
0^~, O2CH, and pyridine ligands, and the total IR density of states. All of the 
projected IR density of states have the same scale as the total IR density of states. 
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Fig. 5. Calculated total and projected Raman vibrational spectra for the same 
monomer used for the IR spectra. All of the projected Raman intensities have the 
same scale as the total Raman intensity. 
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